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We present results on the scalar strangeness and charm contents of the nucleon and of the dis- 
connected contributions to the nucleon spin. These are obtained on np = 2 non-perturbatively 
improved Sheikholeslami-Wilson configurations at a pseudoscalar mass of 290 MeV. We quote 
fr s = m s (N\ss\N) /m^ = 0.070(22) as our preliminary value for the strange quark fraction of the 
nucleon mass and an MS scheme value As = —0.015(10), with as yet unknown systematics, for 
the strangeness contribution to the spin. 
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1. Introduction 

We calculate disconnected scalar matrix elements (N\qq\N) as well as contributions to the spin 
of the nucleon Aq. Of particular interest are currents containing the strange quark q = s, where only 
quark-line disconnected diagrams contribute but we also consider heavier and lighter flavours. 

The scalar matrix element determines the coupling of the Standard Model Higgs boson to the 
quarks within the proton. This then might couple to heavy particles that could be discovered in 
future LHC experiments, some of which are dark matter candidates The fraction of the proton 
mass hin carried by the quark flavour q, fr = m q (N\qq\N) /m^, is particularly important since the 
combination Y, q fr q will appear quadratically in this cross section. It is evident that fj q — ► for 
m q — > and fr oc (N\GG\N)/m^ = fr G for m q — ► oo. We calculate the scalar matrix element for 
quark masses up to the charm quark to confirm this limiting behaviour. It is also not immediately 
obvious whether the charm quark should be considered as a sea quark or not. We remark that in 
general there will be mixing between the dimension four operators GG and m q qq. This deserves 
future study. 

Disconnected contributions to the nucleon structure are also needed with respect to precision 
measurements of Standard Model parameters in pp collisions at the LHC where for instance the 
resolution of a (hypothetical) mass difference between the W + and W~ bosons is limited by the 
theoretical knowledge of the asymmetries between up and down as well as between strange and 
charm sea quark contents of the proton 

The spin of the nucleon can be factorized into a quark spin contribution AE, a quark angular 
momentum contribution L q and a gluonic contribution (spin and angular momentum) AG: 

^ = -AZ+L q + AG. (1.1) 

In the naive non-relativistic SU(6) quark model, AL = 1, with vanishing angular momentum and 
gluon contributions. In this case there will also be no strangeness contribution As in the factoriza- 
tion, 

A£ = Ad + Au + As + --- , (1.2) 

where in our notation Aq contains both, the spin of the quarks q and of the antiquarks q. Exper- 
imentally, As is obtained by integrating the strangeness contribution to the spin structure function 
gi over the momentum fraction x. The integral over the range in which data exists (jc > 0.004) usu- 
ally agrees with zero. For instance a recent Hermes measurement in the region x > 0.02 yields [||] 
As = 0.037(19) (27). This means that non-zero results rely on extrapolations into the experimen- 
tally un-probed region of very small x and are model dependent [Q, Hi. The standard Hermes 
analysis [§] yields As = -0.085(13)(8)(9) in the MS scheme. 

We reported first results and developed the necessary methods in refs. [f7]-§]. These algorith- 

9 

mic studies were performed on rooted staggered rip = 2+ 1 configurations and indicated a value, 
As > —0.02. However, there are unresolved theoretical issues with this fermion approach [ffO], 
so that for the physics study, of which we present preliminary results here, we employ improved 
Wilson fermions that have a meaningful continuum limit. 
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2. Simulation details 

Simulations were performed on 24 3 x 48 and 32 3 x 64 QCDSF configurations of tip = 2 non- 
perturbatively improved clover Wilson fermions with Wilson gauge action at j8 = 5.29 and fc sea = 
0.13632. These values correspond to a pseudoscalar mass mps ~ 290 MeV and a lattice spacing 
a -1 at 2.59 GeV [|Pl|]. The scale was set from the value r ( y 1 at 422 MeV, obtained by chirally 
extrapolating the combination m^ro to the physical point. The spatial lattice extents correspond to 
La rj 1.83fm and La ss 2.43 fm, with a larger 40 3 x 64 volume (La « 3.04 fm « 4.5»ips) still being 



analyzed. We use a modified version of the Chroma code [12]. 

We vary the quark mass parameter Kj oop of the current insertion as well as that of the nu- 
cleon valence quarks JC va i- In particular we combine the values fc = fc sea = 0.13632, K = 0.13609 
and K = fCstnmge = 0.13550, corresponding to pseudoscalar masses «i PS « 290, 440 and 690 MeV, 
respectively. Additional heavier masses are used for the loop. 

Aq and (N\qq\N) are extracted from the ratios of three-point functions to two-point func- 
tions (at zero momentum): 

Re (i^c^(ft)E 1 Tr(iJf- 1 (x,f;x,0ri«»p)) 

For the scalar matrix element we use, r2 pt = r unpo i := (1 + %)/2 and ri oop = 1. For Aq we calculate 
the difference between two polarizations: r2 pt = YjYs^nnpol an d ri oop = Jjjs, where we average 
over all three possible j-orientations. The spin projection operators along the 7-axis read, Pt\ = 
5 (1 ± iyy75)' so tnat m tn i s case > I"2pt = — i(^T — ^ > |)r'unpoi> where we have traded a factor i against 



taking the real part, rather than the imaginary part, of the nominator in eq. (|2. 1|). The variance of the 
above expression is reduced by explicitly using the fact that ImTr(M _1 l) = ImTr {M Jjjs) = 0. 

C a P (tf) denotes the two-point function of the zero momentum projected proton with sink and 
source spinor indices a and /3 and positions tf and t\ = 0. In the limit of large times, tf 3> t 2> 0, in 
the axial case, 

R dis (t,t f )+R con (t,tf)^Aq, (2.2) 

where we have not computed the connected contribution 7? con . This vanishes for strangeness and 
charm contents. In the scalar case the vacuum contribution (0|^|0) = — ^ x Re(Tr(M _1 (x,f;x,f)) 
needs to be subtracted from eq. ( |2.1[ ). We employ sink and source smeared two point functions such 
that this asymptotic limit is effectively reached for t = 4a « 0.3 fm and tf = 8a « 0.6 lfm. We vary 
tf to check this assumption and compute the final results from the tf > 8a data. 

The clover Wilson operator M can be written as 2kM = 1 — kD. We estimate Tr[M _1 T] 
stochastically at the cost of less than 100 preconditioned CG solves on each configuration. The 
noise is reduced by calculating the first two terms of the hopping parameter expansion explicitly, 
which corresponds to multiplying the estimates by (kD) 2 . For details, see [^] and references 
therein. 

For the calculation of Aq we employ the truncated solver method ^] where for the strange 
mass on the 32 3 x 64 lattice, A^i = 730, N2 = 50, n t = 40. For the truncated solves we use the 
smoothly converging even-odd preconditioned CG algorithm while for the runs to full convergence 
we employ the faster BiCGstab2 algorithm. For the scalar matrix element we only use these 50 



3 



Strangeness and charm content of the nucleon 



Gunnar Bali 



A 

z 



z 

v 



Figure 1: The scalar matrix element at Kj oop = 



^strange 



for 8 different sources and the av- 



eraged result. 




Figure 2: The partial sum eq. (2.3) for the scalar 
matrix element at K"i oop = fc va i = ^strange for L = 
32 (red symbols) and L = 24 (black symbols). 



solutions since this is dominated by the gauge noise, not by the stochastic noise. We deviate 
from ref. [^] by seeding the stochastic sources on 8 timeslices rather than on one. This does 
not cause any computational overhead. We then compute the standard point-to-all propagators 
that are needed for the nucleon two point functions for 4 time-separated source points on each 
configuration (requiring 48 solves rather than 12). Exploiting forward and backward propagation 
in time (replacing r2 pt = (1 + /4)/2 by (1 — 74) /2), this then gives us 8 two point (and three point) 
functions per configuration. Indeed, this averaging reduces the errors of As and (,/V|.w'|./V) by factors 
~ 1 / V& which is the maximal possible gain. We display this for the latter example in figure [l| 

In the ongoing analysis of the 40 3 x 64 volumes we will compute the lowest eigenmodes of the 
Hermitian Dirac operator 75M, to further reduce the stochastic noise of the disconnected loops [g], 
to precondition the solver and for low mode averaging of the nucleon two point functions Jl3||. 

In the present study of the L = 24 and L = 32 lattices, where no zero momentum projection is 



performed at the source, we find it worthwhile to rearrange the nominator within eq. (2.1 ): 



I (^C 2pt (y,tf;0,0)TT[M- 1 (x,t;x,t)r])= £ /(x) ~ 4««" mw + • 



(2.3) 



X <-V n - 



x \<x m 



At large x the f(x) ~ e~ mx values will eventually not contribute to the signal anymore but just 
increase the statistical noise. We display the partial sums for the 1.8 fm (L = 24) and 2.4 fm (L = 32) 
lattices for the scalar matrix element at fc va i = fQ oop = fc strange in figure ||. Indeed, at small x max we 
see the expected ^ ax volume scaling. This flattens somewhat around x max « 8a but only saturates 



when the boundaries of the box are hit (x n 



12a and x„ 



16a, respectively). Beyond these 



distances only the lattice "corners" are summed up. This non-saturation means that partial sums 
can only become a permissible method of reducing the noise at much larger spatial volumes. It 
also indicates that finite size effects might still be substantial for the 2.4 fm data. The partial sums 
appear to saturate even more slowly for the case of As. 



3. Results 

In figure ||| we display our La w 1.8fm results on fj for two different nucleon valence quark 
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Figure 3: The proton mass fractions fj for Figure 4: The strangeness fj s for different va- 



quarks q of different masses. 



lence quark masses and volumes. 



masses K = fc s trange (red symbols) and K = 0.13609 (blue symbols), as functions of the quark mass 
m q °c mp S in the loop. The scheme independent combination m q (N\qq\N) encodes the contribution 
of the mass m q of quark q to the nucleon mass. This suggests the normalization with respect 
to the nucleon mass m^. We note that differences between the PCAC and the bare quark masses 
m q = (kt x — JCjL)/ (2a), where K c va \ is obtained from an extrapolation of mp S to zero, only exceed 
the percent level for m q < m s . The right-most data correspond to the charm, the data near nip S ps 
0.48 GeV 2 to the strange quark, fj approaches zero as m q — > and saturates at fj m 1/3 for 
m q > (2/3)m c . Obviously the charm, bottom or top quarks cannot each be made responsible for 
one third of the proton's mass. Moreover, in the heavy quark limit, (N\qq\N) ~ (N\GG\N)/m q , so 
that mixing with the gluonic matrix element needs to be considered carefully. 

In figure |] we focus on the volume and quark mass dependence of fj s . The chiral behaviour, 
varying the valence pseudoscalar mass from 690 MeV down to the 290 MeV that correspond to the 
sea, is well fitted by a constant. So, hopefully, the dependence on the sea quark mass will be weak 
as well. fr s appears to increase on the larger volume, see also figure but to exclude this to be 
just a statistical fluctuation, we will have to wait for the L = 40 analysis. The value obtained on the 
larger volume for the lightest nucleon mass reads, 



fr, 



m s {N\ss\N) 



0.070 ±0.022. 



m N 



Other recent direct calculations resulted in the values fj s = 0.34(5) for anisotropic tip = 2 Wilson 
fermions [14] and fj s = 0.015(28) for np = 2 overlap fermions, fixed to zero topology [13]. An 
indirect determination with rooted staggered fermions, combining chiral condensate and nucleon 
two point function data with the Hellmann-Feynman theorem, suggested fj s = 0.063(11) Jl5|]. 

The valence quark mass and volume dependence of As is displayed in figure |5[ Both appear 
to be mild, with the tendency of a bigger —As on the larger volume, in particular for the lightest 
nucleon mass. We quote the value, 



As 



latt . 



-0.020 ±0.013. 
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Figure 5: Valence quark mass dependence of As latt on the La = 24a 1 .8 fm and La = 24a w 2.4fm lattices. 



that we obtain on the large volume for the lightest nucleon mass. Note that this number applies to 
the lattice scheme and needs to be multiplied by a renormalization factor that we expect to be close 
to 0.76 for a conversion into the MS scheme. Another recent study yielded As = —0.0064(24) [jjjp, 
employing np = 2 anisotropic Wilson fermions. We remark that our value at K" va j 
— 0.0187(54) and, therefore, differs from zero by 3.5 standard deviations. 



^strange reads 



4. Outlook 

At present, the La « 3 fm 40 3 x 64 volumes are being analyzed. Disconnected contributions 



to form factors []16j at non-vanishing momentum transfer and to moments of partem distribution 
functions are also of big phenomenological interest and will obviously extend the present study. 
The long term goal of this project is to go to large volumes at the physical sea quark mass. 
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